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QQ I Abstract 

(N ' 

r^ ■ Quantum mechanics requires the operation of quantum computers to be unitary, and 

(^ ! thus makes it important to have general techniques for devefoping fast quantum al- 

t^^ ' gorithms for computing unitary transforms. A quantum routine for computing a 

generahzed Kronecker product is given. Apphcations include re-development of the 



networks for computing the Walsh-Hadamard and the quantum Fourier transform. 
New networks for two wavelet transforms are given. Quantum computation of Fourier 
P5 ■ transforms for non-Abelian groups is defined. A slightly relaxed definition is shown to 

^ . simplify the analysis and the networks that computes the transforms. Efficient net- 

qh| works for computing such transforms for a class of metacyclic groups are introduced. 

A novel network for computing a Fourier transform for a group used in quantum error- 
correction is also given. 



1 Introduction 

The quantum computational version of the discrete Fourier transform is without doubt the 
most important transform developed for quantum computing so far. It is in the heart of all 
quantum computational issues discussed until now. All main quantum algorithms, including 



Shor's celebrated factoring algorithm |29] and Grover's searching algorithm [21], use it as 
a subroutine. All known relativized separation results for quantum computation are based 
on quantum algorithms that use the discrete Fourier transform p, H, [3^, Q. Fundamental 
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concepts of quantum error-correction rely on it; see for example [^ ^, |T^, |HT], ^. It is, 
in conclusion, the most important single routine for obtaining an insight in the previous 
work done in quantum computing. Still, this seemingly simple transform is not yet fully 
understood. 

When referring to the discrete Fourier transform, one often does not refer to a single trans- 
form, but rather to a family of transforms. For any positive integer n and any ra-dimensional 
complex vector space Vn-, one defines a discrete Fourier transform F^ (see for example 



More generally, given r positive integers {nj}[^]^, and r complex vector spaces {V^}[^]^, where 
Vi is of dimension rij, one defines a discrete Fourier transform, denoted F„^ ® ■ ■ ■ ® Fn^, for 
the tensor product space V^i (S) ■ ■ ■ (S> K- (See Section |^ for details.) 

Earlier, the quantum versions of the discrete Fourier transforms were defined, but efficient 
quantum networks where only known for few of them. Currently, efficient networks imple- 



menting Fn exactly have been found for all smooth integers n |Tl|, |1^, |I6| ^, |2^, where n 
is considered smooth if all its prime factors are less than log'^(n) for some constant c [29|. 
Furthermore, any discrete Fourier transform can be efficiently approximated to any degree 
of accuracy by some quantum circuit |0, Q . 

Common for the networks discussed above are that their description has been taken from 
the point of view that the transforms were to be implemented by quantum networks. In this 
paper, we discuss quantum computation not as opposed to classical (perhaps parallelized) 
computation, but more as a variant. We believe that the fundamental object is the unitary 
transform, which then can be considered a quantum or a classical (reversible) algorithm. 
With this point of view, the problem of finding an efficient algorithm implementing a given 
unitary transform U, reduces to the problem of factorizing U into a small number of "sparse" 
unitary transforms such that those sparse transforms should be known to be efficiently 
implementable. As an example of this, we show that the quantum networks implementing 
the quantum versions of the discrete Fourier transforms can be very easily derived from the 
mathematical descriptions of their classical counterparts. 

More generally, we consider a new tool for finding quantum networks implementing any given 
unitary transform U. We show that if U can be expressed as a certain generalized Kronecker 
product (defined below) then, given efficient quantum networks implementing each factor in 
this expression, we also have an efficient quantum network implementing U. The expressive 
power of the generalized Kronecker product includes several new transforms. Among these 
are the two wavelet transforms: the Haar transform [^ and Daubechies' V^ transform |H|, 
and we are thus able to devise new quantum networks that compute these transforms. 

There exists a group theoretical interpretation of the discrete Fourier transforms which es- 
tablishes a bijective correspondence between this family of transforms and the set of finite 
Abelian groups. To further demonstrate the power of the generalized Kronecker products, 
we use them to give a simple re- development of the quantum Fourier transforms for Abelian 



groups. More interestingly, the Fourier transforms can be generalized to arbitrary finite 
non-Abelian groups (see for example ||25[ for an introduction), and we give a definition of 
what it means that a quantum network computes such transforms. Moreover, we give a 
slightly relaxed definition where we only compute a Fourier transform up to phase factors. 
Classically, the idea of relating a Fourier transform for a group to one of its subgroups has 
proven to be very useful 



25 



and we show that this carries over to quantum computers. 
We apply these ideas to give new networks for quantum computing Fourier transforms for 
the quaternionic group and for a class of metacyclic groups. 

Since Shor demonstrated that quantum error-correction is possible by given an explicit nine- 
bit code [0, several new classes of quantum codes have been developed. See for example 

Many of these are stabilizer codes, which 



10, 18, 31, 32| for some of the many results. 



are subgroups of a certain non-Abelian group En (defined in Section |^). For that group, we 
also give a simple and efficient network for computing a Fourier transform, again using the 
framework of generalized Kronecker products. 

Independently of this work, Robert Beals has found quantum networks implementing Fourier 
transforms for the symmetric groups [0]. A challenging open question related to that result 
is whether it can be used to find a polynomial quantum circuit solving the famous graph 
isomorphism problem. 



2 Generalized Kronecker products 

All matrices throughout this paper are finite. Matrices are denoted by bold capital letters 
and tuples of matrices by calligraphic letters. Indices of tuples and row and column indices of 
matrices and vectors are numbered starting from zero; the {i,j)-th element of A is referred 
to as aij. A single integer as subscript on a unitary matrix denotes its dimension, e.g., we 
let Iq denote the [q x q) identity matrix. The transpose of A is denoted by A*. Recall that 
a square matrix is unitary if it is invertible and its inverse is the complex conjugate of its 
transpose. The complex conjugate of a number c is denoted by c. 



Definition 1 Let A be a (p x g) matrix and C a (A; x /) matrix. The left and right Kro- 
necker product of A and C are the [pk x ql) matrices 



Acoo Acoi . . . Aco,/-i 

Acio Acii ... Aci,/_i 

_Acfc_i,o Acfc-1,1 . . . Ack-ij-i 
respectively. 



and 



CbiQ\^ Cbiiy^ 

.tip-i,oC ap_i^iC 



0'p-l,q-lC 



n 



We denote the left Kronecker product by A ®i C and the right Kronecker product by 
A ®pi C. When some property holds for both definitions, we use A ® C. Note that the 
Kronecker product is a binary matrix operator as opposed to the tensor product which is 
binary operator defined for algebraic stuctures like modules. The Kronecker product can be 
generalized in different ways; see for example [|l^] and |2^. In this paper, we use (an even 
further generalized version of) the generalized Kronecker product discussed in [jl^, defined 
as follows. 



Definition 2 Given two tuples of matrices, a fc-tuple A = {A^)^Iq of (p x q) matrices 
and a g-tuple C = {C^)1Zq of [k x /) matrices, the generalized right Kronecker product is 
the {pk X ql) matrix D = ^ 0^^ C where 



"■ij (^uk+v,xl+y ^ux'^vy 



with 0<u<p, 0<f<A;, 0<a;<g, and < y < I. 
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The generalized right Kronecker product can be found from the standard right Kronecker 
product by, for each sub-matrix a^xC in Definition |I] substituting it with the following sub- 
matrix 






'^ux'^01 
^ux^ll 



^ux ^k-l,0 ^ 



k—l„x 



■'k-lA 



'^ux^O,l-l 
^ux^ll-l 



"'ux '^k-lJ-1. 



The generalized left Kronecker product is the {pk x ql) matrix D = ^ ^l C where the 
(i, j)-th entry holds the value 



lAlri n Ct" 



up+v,xq+y 



^vy^ux 



with 0<'U<A;, 0<f<]9, 0<x</, and < y < q. 

As for standard Kronecker products, we let A^C denote either of the two definitions. 
If the matrices A* = A are all identical, and also C* = C, the generalized Kronecker product 
A<^C reduces to the standard Kronecker product A ® C. Denote by ^ ® C the generalized 
Kronecker product of a fc-tuple ^ of (j9 x q) matrices, and a g-tuple C of identical {k x /) 
matrices C. Denote A ® C similarly. 
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To analyze generalized Kronecker products we need the shuffle permutation matrix of di- 
mension {mn X mn), denoted Hmn as shorthand for Il(^rn,n), defined by 



TTr. 



''d7i+e,c('rra+e' 



^de'Sd'e, 



where < d,e' < m; < d',e < n, and 6xy denotes the Kronecker delta function which is 
zero ii X ^ y, and one otherwise. It is unitary and satisfies n~^ = n^„ = Unm- 

Given two tuples of matrices, fc-tuple A = (A*)fro^ of (p x r) matrices and fc-tuple C = 
(C*)fro^ of (r X q) matrices, let AC denote the fc-tuple where the i-th entry is the {p x q) 
matrix A*C*, < i < k. For any k-tnple A of matrices, let Diag(^) denote the direct 
sum ^i^Q A* of the matrices A*^, . . . , A'^"^. The generalized Kronecker products satisfy the 
following important Diagonalization Theorem of [1^ . 



Theorem 4 [Diagonalization Theorem] Let A = (A*)j^q be a /c-tuple of {p x q) matrices 



and C 



[C*)f=o a g-tuple of {k x /) matrices. Then 



A®rC= [Upk Diag(^) Ukq) X Diag(C) 
A®lC = Diag(^) X (nfc,Diag(C) n^. 



(1) 
(2) 



Corollary 5 Let A = (A*)j^q be a A;-tuple of {p x q) matrices and C = {C)1^q a g-tuple 
of {k X /) matrices. Then 

A®rC = Upk[A^LC)Ulq 
A(»lC = Ukp(A®RC)Uql. 



Until now, we have not assumed anything about the dimension of the involved matrices. 
In the next theorem, we assume that the matrices involved are square matrices. The theorem 
is easily proven from the Diagonalization Theorem. For any A;-tuple A = {A^)^~q of invertible 
matrices, let A~^ denote the fc-tuple where the i-th entry equals the inverse of A*, < i < k. 

Corollary 6 Let A, C be ?7i-tuples of {n x n) matrices, and V, S be ra-tuples of {m x m) 
matrices. Then 

(3) 



AC®V£ = i^A® Irn) X [C0Vj X [I.^ ® £j . 
Furthermore, if the matrices in the tuples A and C are invertible, then 

[a ®r Cy^ = U„,m {C-^ ®fl A-^) n„„ = c-^ ®L A~^ 

[a ®L C) "' = Umn {C-^ ®L A-^) n„^ = C-^ 0R A-^ 

If A and C are unitary, then so is ^ C. 



3 Quantum routines 

In this section, we give a method for constructing a quantum network for computing any 
given generalized Kronecker product. A primary application of this method is as a tool 
to find a quantum network of a given unitary matrix. As two examples, we use it to de- 
velop quantum networks for computing two wavelet transforms, the Haar transform and 
Daubechies' D^ transform. 

As our quantum computing model, we adopt the now widely used quantum gates arrays [|l], H. 
Let T : \u, v) \-^ \u,v Q) u) denote the two-bit exclusive-or operation, and U the set of all one- 
bit unitary operations. Following 0] , by a basic operation we mean either a U operation or 
the r operation. The collection of basic operations is universal for quantum networks in the 
sense that any finite quantum network can be approximated with arbitrary precision by a 
quantum network Q consisting only of gates implementing such operations [13, 

Define the one-bit unitary operations 
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Given a unitary matrix C, let A((j, x), {k, C)) denote the transform where we apply C on the 
k-th register if and only if the j-th register equals x. Given an n-tuple C = (C*)"rQ^ of unitary 
matrices, let A{{j,i),{k,C% denote A((j>- 1), (/c, C"-^)) ■ ■ ■ A((j, 0), (fc, C°)). Given a 
/c-th root of unity, say cu, let $(0;) denote the unitary transform given by \u)\v) i— > uj^'"\u)\v). 
If the first register holds a value from Z„, and the second holds a value from Z^, then 



$r 



'{n,m)i^) can be implemented in 9([logn] [logm]) basic operations ||TT], |T2|, ^] . 

Quantum shuffle transform. For every m > 1, let the operation DM^ perform the 
unitary transform \k)\0) ^-^ \k div m)\k mod m). Let SWAP denote the unitary transform 
\u)\v) ^— *> |f)|M). Then Hmn can be implemented on a quantum computer by one apphcation 
of DMm, one swap operation, and one application of DM~^, 

n^„ = dm;1 swap dm^. 

Quantum direct sum. Let C be an n-tuple of {m x m) unitary matrices. It is not difficult 
to see that Diag(C) can be implemented as follows, 

Diag(C)=DM-i A((1,2),(2,C*)), DM^. 

In general, the time to compute the direct sum is proportional to the sum of the computation 
times of each of the conditional C* transforms. However, if parts of these can be applied in 
quantum parallel, this improves the running time. 

Quantum Kronecker product. Let A be an m-tuple of {n x n) unitary matrices and C 
an ra-tuple of {m x m) unitary matrices. By the Diagonalization Theorem, the generalized 
Kronecker product can be applied by applying two direct sums and two shuffle transforms. 
Removing cancelling terms we get 



A0rC = BM-^ A((2, t), (1, A^)), A((l, z), (2, C^)), DM^ 
A^lC = DM;1 A((1, ^), (2, A^)), A((2, 2), (1, C% DM,. 



(4) 
(5) 



Thus, an application of a generalized right Kronecker product can be divided up into the 
following four steps: in the first step, we apply DM^- In the second step, we apply the 
controlled C* transforms on the second register, and in the third step, the controlled A* 
transforms on the first register. Finally, in the last step, we apply DM^^ to the result. 

Example 7 Let ^ be a 4-tuple of (2 x 2) unitary matrices, and C a 2-tuple of (4 x 4) unitary 
matrices. The generalized Kronecker product A ®ij C can be implemented by a quantum 
network as follows. 



MSB 
LSB 



^rrt 



a^Ha 



^&^ 



A' 



A= 



■MSB 
■LSB 



The dots and the circles represent control-bits: if the values with the dots are one, and if the 
values with the circles are zero, the transform is applied, otherwise it is the identity map. 
The least (most) significant bit is denoted by LSB (MSB). Note that most of the transforms 
are orthogonal, and thus the gates commute. Note also that, following the ideas of Griffiths 
and Niu in |l^ , a semi-classical generalized Kronecker product transform can be defined. D 



3.1 Quantum wavelet transforms 

A main application of the Diagonalization Theorem is as a tool to find quantum networks 
computing large unitary matrices. Suppose we have factorized a unitary matrix U via a 
generalized Kronecker product into a product of some simpler matrices. Then, if we have 
quantum networks for computing these simpler transforms, we also have a quantum network 
for computing U by applying the methods developed above. We give two examples of this 
technique; in both cases implementing a wavelet transform. 



Example 8 The Haar wavelet transform |]2^, H, can be defined using the generalized Kro 
necker product as follows. 

H2 = W 

Hgn+l = 112,2- X ((H2n,l2n) 0/? W) , ?2 = 1 , 2 , . . . 



(6) 



Applying the decomposition of the generalized right Kronecker product given in Equation (4), 
we immediately obtain an efficient quantum circuit for computing the Haar transform. 
Let 82"+! denote the bit-shift transform given by |6„ . . . 60) '-^ \bobn . . . &i). This transform 
efficiently implements 112,2" • For n = 3, the quantum circuit defined by Equation (|) is given 
below. Here the two S transforms are the bit-shift transforms of the appropriate dimensions. 
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w 



s - 



n 

By Equation (|^), we can rewrite the recursive definition in Equation (P) as 

H2n+i = H2,2" X ((H2n,l2") ®Rh) X (I2" (8>ij w) . 

We refer to the right-most factor in this factorization as the scaling matrix of dimension 
(2"+^ X 2"+^) for the Haar wavelet transform. In general, given any family of unitary matri- 
ces {D2i}i>io, define a family of unitary transforms {U2'}i>io as follows, 

LI 2*0 = ■'-'2'0 
U2"+»o ^ J-J-2 2"+'0-i X {y\J2n+io-l,i-2n+iQ-lj'S)Ri-2) ^ ^2"+'0i n ^ if Zy . . . 



We refer to the Dj matrices as a family of scaling matrices, and the family of Uj matrices 
as a wavelet transform. Suppose that we have a family of efficient quantum networks for 
computing a given family of scaling matrices. Then, as in Example ^, we also have efficient 
quantum networks for computing the associated wavelet transform. The next example gives 
a factorization of the scaling matrices used in Daubechies' D^ wavelet transform fl^ 

Example 9 Let m > 4 be an even integer, and let 

, 3±V3 , , 1tV3 

Daubechies' D^ scaling matrix [|14| of dimension {m x m) is the matrix with 

^ kj_i^x i is even 

dij — 



where 



X 



-iyk2+i-j-x i is odd, 



4 a i > m — 2 and j < 2 
otherwise. 



and fc; = if / < or if / > 3. 

Let Pm be the {m x m) permutation matrix which subtracts two if the input is odd, i.e., 
Pij = 1 a i = j and i is even, or if i + 2 = j (mod m) and i is odd. Let Co and Ci denote 
the two one-bit unitary operations, 

Co = 2(1' ~lA and C, = 1 N^^ ^]A. 

The scaling matrix D^ can then be factorized using two Kronecker products 

D^ = (l„/2 ^R Ci) X P„ X (l„/2 ®R Co) . (7) 

Set n = [logTTi]. The permutation transform P^ can be implemented in Q{n) basic op- 
erations |34|. Each of the other two factors on the right hand side of Equation (|^) can be 



implemented in one basic operation. Thus, D^ can be implemented in B(n) basic opera- 
tions. We remark that we have not been able to find this factorization of D^ elsewhere in 
the literature — despite the fact that using it compared to straightforward use of D^ saves 
m additions in the classical case. Furthermore, note that flm/2 ®r Cij x (1^/2 ®r Co) = 
Im/2 ®i? W, which is the scaling matrix used in the Haar transform. D 
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4 Group representations and quantum Fourier 
transforms 

In the rest of this paper, G will denote a finite group, written multiplicative with identity e, 
and T} the order of G. Let CG denote the complex group algebra of G. Let ;Btime denote the 
standard basis of CG, that is, {gi, . . . , Qr,}, and let (m, v) = SgeG u{g)v{g) denote the natural 
inner product in CG. Let GLrf(C) denote the multiplicative group of {d x d) invertible 
matrices with complex entries. We start by reviewing some basic facts from the theory of 
linear representations of finite groups. For a general introduction to group representation 



theory, see for example [|T3| or ||2l 

A complex matrix representation p of G is a group-homomorphism p : G ^ GLd(C). The 
dimension d = dp is called the degree or dimension of the representation p. Two representa- 
tions, pi and p2, of degree d are equivalent if there exists an invertible matrix A G GLrf(C) 
such that p2{g) = A~^pi{g)A for all g E G. A representation p : G ^ GLd(C) is irreducible 
if there is no non-trivial subspace of C^ which is invariant under p{g) for all g E G, and it 
is unitary if p{g) is unitary for all g E G. For every representation there exists an equiva- 
lent unitary representation. Up to equivalence, there are only a finite number of irreducible 
representations, say u, of G. This number equals the number of distinct conjugate classes 
of G. 

Let TZ = {p^, . . . , p^} be a complete set of inequivalent, irreducible and unitary represen- 
tations of G with di equal to the degree of p*. For any representation p E TZ, the vector 
Pki G CG defined by considering the [k, /)-th entry of p{g) for each g E G is called a matrix 
coefficient of 7^. The inner product of two matrix coefficients of 7^ is non-zero if and only if 
they are equal. For each matrix coefficient pki, let bp^k,i denote the normalized matrix coeffi- 
cient, and let Bf^^q = {bp,k,i} denote the set of orthonormalized matrix coefficients. Since one 
can show that the degrees di of the representations pi ElZ satisfy the relation Yh=i df = rj, 
it follows that i3frcq is an orthonormal basis of the vector space CG. 

The linear operator Fq on CG which maps a vector v G CG given in the standard basis -Btime 
to its representation v G CG in basis iSfreq is called the Fourier transform for CG on TZ. Each 
entry of -0, denoted v{pki) or just pki, is called a Fourier coefficient of v (on 7Z). 

In the recent years, many new exciting results have been found concerning the computation of 
Fourier transforms for finite groups on classical computers — see p5[ for a nice survey. In this 



paper, we consider the computation of Fourier transforms on quantum computers. Since 
quantum mechanics requires the operation of the computer to be unitary, our definition of a 
Fourier transform given above is slightly more strict than the most common used definitions 
for the classical case. 
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We now define what it means that a quantum circuit computes a Fourier transform. Let Fq 
be a Fourier transform for CG on IZ. Let -Etimc '■ -Btimo —^ ^t) and -Efreq : -Bfreq — ^ ^»7 be 
two bijections. These functions induce an ordering on ;Btime and ;Bfreq, respectively. Let E : 
;Btime U ;Bfreq ^ Z^ deuote the extension of -Etime and -Efreq- We say that E is an encoding for 
the hnear transform Fq- With respect to E^ Fq can be viewed as a matrix F^ in GL^(C). 
This matrix is unitary by construction, and thus there exists a quantum circuit computing 



it ||15[. We say that the circuit computes Fq with respect to E. 

Given a fc-tuple of complex numbers of unit norm, (0i)^=i, let = diag((/)j) G GLfc(C) denote 
the unitary diagonal matrix with (pi at the i-th diagonal entry. Let Fq be a Fourier transform 
for CG on TZ, E an encoding for Fq, and F^ the resulting unitary matrix. We say that a 
quantum circuit computes Fq up to phase factors (with respect to E) if there exists a unitary 
diagonal matrix G GL^(C) such that the circuit computes F^ = 0Fg. Given a network 
that computes Eg up to phase factors, we can obtain a quantum circuit for computing Fq 
exactly by applying first F^ and then the unitary transform 0~^ = 0*. 

Note the dependencies of the set TZ and the encoding E of Stimc and Sfreq in the above 
definitions. A Fourier transform for CG is defined only with respect to TZ. A quantum circuit 
computing the Fourier transform is, in addition, defined with respect to an encoding E of 
the basis-elements in Stime and i3freq- 

The quantum computation time of Fq with respect to TZ and the encoding E is defined as 
the minimum number of basic operations in any quantum circuit computing F^, and it is 
denoted by QT(G)(7?., £"). The quantum computation time of a Fourier transform for CG, 
denoted QT(G), is defined as the minimum of QT(G)(7^, i?) over all possible choices of IZ 
and E. 

In the rest of this paper, TZ denotes a complete set of inequivalent, irreducible and unitary 
representations of G. 



5 Quantum Fourier transforms for cyclic groups 

As an introductionary example, consider the problem of quantum computing the discrete 
Fourier transform. We start by developing an efficient quantum routine for computing the 
discrete Fourier transform using the generalized Kronecker product discussed in Section ^ 
Then, we review a group theoretical interpretation of the transform which relates it to the 
cyclic groups. 
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The discrete Fourier transform for a quantum computer is defined as follows. For any positive 
integer n, let 

1 n—l 
Fnk) = -=E^n1l/), (8) 

for each x = 0, . . . ,n — l, where Un = exp{2TT\^—l/n) is the principal n-th root of unity. The 
unitary Fourier transform Fnm can be defined from F^ and F^ using a generalized Kronecker 
product 

Fnm = n„„ X (F„ (g)L Im) X ((D^„)^Jo^ (g)L I-m) X (l„ ®L F„) (9) 

where D^^ = diag(c<j'*') for < s < m and u = Unm- Equation (^ is referred to as a radix^n 
splitting in [Q, where a proof of the identity can be found. 

Equation (^) gives an efficient quantum routine for computing Fnm from F„, F,„, and D^^. 
Interestingly, the resulting routine obtained this way is the same as the one found by 
Cleve |]TT| using a direct method. The transform ((D^^)^ ®l I^) is a special application 



of the $ transform defined in Section ^ and is thus easily applied. For powers of 2, the 



computation of F2n uses 6(n^) basic operations [12|. 

We now review the well-known group theoretical correspondence to the discrete Fourier 
transform. Let G = Z„ be the cyclic group of order n. For Abelian groups all irreducible 
representations are one-dimensional, and hence equivalent representations are equal. There 
are n distinct representations, IZ = {C°, . . . , C""^}; given by 

C{j) = K] for every j G Z„. 

The collection of normalized matrix coefficients are i3freq = {^(o 1 1, . . . , 6^n-i 1 1} where 
(V,i,i>i) = -^^n for all b^^^l^l G Sfreq and j G Stimc- Hence, for all b^^^l^l G Bfreq we 
have 6^«^i_i = -^ EjeBtime '^n 3- Thus, by choosing the encoding E given by -Etime(j) = j and 
-Efreq(&C%i.i) ~ ^^ respectively, the quantum circuit defined by Equation (H) is seen to compute 
the Fourier transform for the cyclic group Z„ with respect to E. We remark that it is possible 
also to give a group theoretical interpretation of the decomposition given by Equation (P); 



see for example [25| for details. 



5.1 Direct product groups 

Suppose we are given quantum networks (as black-boxes) for quantum computing Fourier 
transforms for the group algebras CGi and CG2- Consider the problem of computing a 
Fourier transform for the direct product group algebra CG = C(Gi x (^2)- Classically, this 
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problem has a very simple solution. In this section, we show that this carries over in the 
quantum circuit model. 

Let Gi and G2 be finite groups of order 771 and 772, respectively. Our first task is to establish 
a specific isomorphism ip between the algebras CGi x CG2 and C{Gi x 6*2) • Let Bl:^^^ denote 
the standard basis of CGj, i = 1,2, and let Btime = {{91,92) '■ 9i ^ Gi, i = 1,2} denote the 
standard basis of C(Gi x ^2)- Let CGi ® CG2 denote the tensor product algebra of CGi 
and CG2, and ip : CGi ® CG2 -^ C(Gi x G2) the natural algebra isomorphism defined by 

<^{9i ® 92) = {91,92) {91 ®92^ B^^^^ O Bl^^ . 
With these definitions, we can write 

Bti^e = ^(Stimc®^Le)- (10) 

Let TZi be a complete set of inequivalent, irreducible and unitary representations of Gi, 
i = 1,2. We need the following lemma from representation theory. 

Lemma 10 Let Gi and TZi be given as above, i = 1,2. Then 

7^ = 7^l (g)^ 7^2 = {p^ (g)^ p2- Pi^ TZi, i = l,2} 
is a complete set of inequivalent, irreducible and unitary representations of G = Gi x G2- 

Let ;B|rgq denote the set of orthonormalized matrix coefficients of TZi, "^ = 1, 2, and B^^eq the 
set of orthonormalized matrix coefficients of 7^, where 7^ is given as in Lemma |TU[ By that 
lemma, it follows by straightforward algebra that 

i3freq = V'(i3Lq®i3|J. (11) 

Having established the above isomorphism, we now state the main result from representation 
theory to be used in this section. Essentially, it reduces the problem of computing a Fourier 
transform for C(Gi x G2) to those of computing Fourier transforms for CGi and CG2. 

Theorem 11 Let Fi and F2 be Fourier transforms for CGi and CG2 on T^i and 7^2, respec- 
tively. Define the linear transform Fq over CGi ® CG2 by 

FG{9i^92) = F^{g^)®F2{92). 
Then Fq = ipFQ{p~^ is a Fourier transform for CG = C{Gi x G2) on 7^ = 7?.i ^r 7^2- 
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This reduction is, however, only as abstract computations over vector spaces. To give a 
concrete quantum circuit for computing the Fourier transform for the product group, we 
also need to consider the choices of bases for the involved transforms in the reduction. 
Let Ei be an encoding for Fi and Fj the corresponding matrix representation of Fi, i = 1,2. 
In ket-notation, the transform F^ reads 

1^1)1^2) I ' (Fi|^i))(F2|^2 

for all c/i ® c/2 e Bl^^^ ® Bl^^. 

Define the bijections Etime : Bume -^ "^^1^2 and Efreq : Bf^eq -^ ^.^^^3 by 

^timc(</'(^l ® ^2)) = V2Ei{gi) + ^2(^2) , . 

Let E denote the extension of -Etime and -Efreq- With respect to the encoding E, Fq as defined 



in Theorem |Tl] has the matrix representation 

Fg = Fi®/jF2. 

Applying Fg is thus done by applying Fi on the most significant bits, and F2 on the least 
significant bits. We have shown 

Theorem 12 Let Gi and G2 be groups of order r/i and 772, respectively. Suppose we have 
quantum networks Fi and F2 for computing Fourier transforms for CGi and CG'2 with re- 
spect to the encodings Ei and E2, respectively. Then the following quantum circuit computes 



a Fourier transform for C(G'i x G2) with respect to the encoding E as defined in Equation ([T^) 



^i^^Fi 
m 



F,- 



As an application of this, consider the Walsh-Hadamard transform p3| , p5|] , defined as follows. 
For any positive integer n, let 



2"-l 



w2n|a;) = 4^^(-i)i:r:'o'-»^»li/). 



On 

^ j/=0 

for each x = 0, . . . , 2" — 1, where x = Xn-i ■ ■ -Xq and y = yn~i ■ ■ -Vo- This unitary transform 
can also be defined using the standard Kronecker product as follows 

W2 = W, W2n+i =W®^W2n, n= 1,2,... (13) 
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Appealing to the generalized Kronecker product routine in Section |^, we immediately obtain 
the well-known method for computing the Walsh-Hadamard transform W2n on a quantum 
computer ||T^: apply the transform W on each of the n qubits. 

It is easy to check that W is the Fourier transform for the cyclic group Z2 of order two. 
Thus, by Theorem IT2|, we have the well-known fact that the Walsh-Hadamard transform 



coincides with the Fourier transform for the Abelian group Zg. This transform has been 
extensively used in quantum algorithms, for example by Deutsch and Jozsa |TB|, Simon [^| 



and Grover [^, §]. One of its advantages is that it can be computed in only Q{n) basic 



operations |jT6| 



Now one might ask if a similar statement holds for subgroups in general. That is, if i7 ^ G is a 
subgroup and we encode elements of G using two registers, the first for coset representatives, 
the second for elements from if, to what extent do the gates then need to involve both 
registers? We consider this question in the next section, and give an answer to it for some 
classes of non- Abelian groups. 

In this section, we have carefully distinguished between isomorphic vector spaces in order to 



prove Theorem |T2|. In the following, we will relax slightly upon this to avoid cumbersome 
notation. Let U and V be any two inner product spaces of dimension m and n, respectively, 
with orthonormalized bases {«!,...,«„,} and {wi,...,f„}, respectively. Then the tensor 
product U ®V and the vector space spanned by {{ui^Vj) : 1 < i < m,l < j < n} are iso- 
morphic under the natural isomorphism Lp given by '^{ui ®Vj) = {ui, Vj). When appropriate, 
we will not distinguish between Ui ® Vj and {ui,Vj) in the rest of this paper. Note that the 
set {ui C?> Vj} is an orthonormalized basis for U (E)V. 



6 Adapted representations 

In the previous section, we related a quantum Fourier transform for G = Gi x 6*2 to quantum 
Fourier transforms for Gi and G2. In the classical case, relating a Fourier transform of a 
group to a Fourier transform to one of its subgroup has shown to be very useful; see for 



example |2^ and the references therein. The main ideas in this approach are factorization 



of the group elements, and the use of an adapted set of representations. For example, in 



Section pTT| , we used the factorization {gi, (72) = (fi'i, 62) ■ (ei, (72)5 where e, denotes the identity 



of Gi, z = 1,2. 

For any subgroup H ^ G and any representation p of G, let p I H denote the representation 
of H obtained by restricting p to H. The representation p | if is unitary but not necessary 
irreducible. Recall that we in this paper assume that all representations are irreducible and 
unitary. 
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Definition 13 Let H ^ G he a subgroup, and 7^ be a complete set of representations of G. 
Then TZ is called H -adapted if there is a complete set 7?.^ of representations of H such that 
the set of restricted representations (TZ I H) = {p I H : p E 71} is a set of matrix direct 
sums of the representations in T^-'^. D 

The set TZ is said to be adapted to a chain of subgroups if it is adapted to each subgroup in 
the chain. Adapted representations always exist. 

Let H ^ G he a subgroup of order m, and T a left transversal for H in G. Let TZ^ he a 
complete set of representations of H, and let TZ he a complete set of representations of G 
that is iJ-adapted relative to 7Z^ . Let B^^^ and Bf^eq denote the collections of normalized 
matrix coefficients for H and G, respectively. Let p G 7^ be a representation of degree d. 
The matrix coefficient pki G CG can be written as a linear sum of the basis-elements Bume- 

Pki = J2 Pkiia) 9 = J2 Y. Y.Pki{'t)pii{h)th = Y.Y.Pki{m Y. Pa{h)th\. 

g&G t€Th£Hi=l t&Ti=l \ h£H ) 

Since IZ is if -adapted by assumption, p is a matrix direct sum of representations in IZ^ . 
Therefore, either puih) = for all h ^ H, or there exist p' G TZ^ of degree d' and 1 < i', I' < d' 
such that pa{h) = p'iniih) for all h E H. In the former case, let p^,;, (and hpi^ii^ii) denote the 
zero vector in 'CH . Then we have 



Pki 



T.T.p^m{T.p^'Ah)th\, 



so 



bpAi = EE^M.wf E p^'Ah)th) = EE\/f V.wf E K'.'Ah)th]. (i4) 



A Fourier transform, F^, is a change of basis in CG from the standard basis to a basis of 
normalized matrix coefficients. Let Fu he the Fourier transform for Cif on IZ^ . For obtain- 
ing an if-adapted method for computing Fg, consider the complex vector space spanned by 
the basis 

T®Bl^, = {t®h : teT.h^Bl,^:^. 

This vector space is clearly isomorphic to CG under the natural map (^ : (T Cg> B^^^^) -^ 
(-Stimc) given by ip(t ® /i) = th. Here, and in the rest of this paper, (■) means span(-). 
Another basis is 

Btemp =T(g) i3ff,q = {t (g) V,^',/' : t^T, bp> ,, ,, G Sffeq}, 
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and using 99, Equation (jT^) reads 



igTi=l ' ^ 



Let l^ : (Sfreq) ^ (-Btcmp) denote the transform 



^ : bp^k,i I — > XI Yl y -7;bp,kA't) t ® V,«',''- 



ter 

By construction of V, 

(/ ® Fh) oip-^ = VoFG 

as illustrated in the following commutative diagram 



/®-Fh 



^G 



(r^i^ff,,) ^— - (i3freq) 



V^ 



Since ^9 is an isomorphism, and F^ and Fq are unitary, ^ is invertible. Let U : (-Btcmp) ^ 
(-Sfrcq) denote the inverse of V, that is, 

^ • Xm \j-Tbp,k,i{t) t ® 6p/,i/,z' I — > bp^k,i (15) 

which maps a vector w G (iStemp) given relative to basis Btemp to its representation v G CG 
relative to basis Sfreq. Hence, we have factorized the Fourier transform Fq into a product of 
three unitary transforms, 

Fg = U o {I Fh) o ip~\ (16) 

A quantum implementation of the adapted method for computing a Fourier transform can 
be obtained as follows. Given a vector v G CG, let Vt G CG denote the vector which is 
non-zero only on the coset tH, on which it is given by Vt{h) = v{th) for all h & H. Initially, 
we hold the superposition v = Z^geBtimo ^(S')!^') ^^^ ^^ want to compute the superposition 
V = J2bi&B[,c '^ibi)\bi). The quantum routine consists of three steps. In the first step, we 
apply ip~^, computing 

v= E v{g)\g)^Y: E v{th)\t)\h) = j:\t)( E v,{h)\h)]. 
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Then, we apply a quantum Fourier transform F^ with respect to T?."^ to the second register, 
producing 



I ircq I rrcq 



Finally, in the third step, we apply the linear transform U given by Equation (|T3|), producing 

E v{h,)\h) = v. 



heBi, 



rcq 



a'c 



(-1: 

1 



The transform U is unitary since the first two steps (that is, (/ Fh) oLp~^) are unitary and 
the composition of the three steps (that is, Fq) is unitary. In the following sections, we apply 
the technique just described to develop quantum Fourier transforms for some non-Abelian 
groups. 

6.1 The quaternionic groups 

The quaternionic group Qn of order An is the group 

Q^ = (^r,c : r^'' = c^ = 1, cr = r^^'-'c, c^ = r''). 

For simplicity, we consider only the case that n is even. The case when n is odd is very 
similar, and in fact slightly simpler. When n is even, Qn has a complete set TZ consisting of 
four one- dimensional and n — 1 two-dimensional representations 

P' = 1 p\r) = p\-c) = 1 ^( x_W 

p\-r) = p\c) = 1 p\-r) = p\-c) = 1 "^ ^""^ ~ [ cJ- 

where 1 <i < n and uj = U2n- 

The group has a cyclic subgroup H generated by r of index two. Let T = {e, c} be a left 
transversal for H in Qn, and write Qn = TH. Let TZ^ denote the complete set of one- 
dimensional representations of H given in Section |^. The set of restricted representations 
of IZ is 

{niH) = {C°, C) U {C' © C'""' : 1 < / < n}, 
so 7Z is if -adapted. Let iStime, -Sfreq, -B/fj^g, and B^^^ be defined as in Section ^. Let 

Stemp = T® Bt^ = {t® 6^.,i,i -.teT, 6^.,i,i G Bt^}. 

A main part of the development of subgroup-adapted Fourier transforms is the determination 
and implementation of the transform U defined by Equation (0). For this purpose, consider 
the matrix coefficient a\^ G CQn- By writing 

teTheH xeZ2n 
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we have that 



Oo-M, 



E 



U'^r^ 



^(e®C) 



Each of the other Fourier coefficients can similarly be written as a linear sum of the basis- 
elements Btemp, 



ba\l,l 


= vie ® C) 


V,i,i 


bcr\l,2 


= (-l)>(c®C'""0 


V,i,i 


ba\2,l 


= vic^C) 


v,l,l 


bcr\2,2 


= vie ® C'"~0 


v,l,l 



^(v9(e®C°)+¥'(c®C°)) 
^(^(e®C°)-¥'(c®C°)) 

^(¥p(e®C)-</'(c®C)) 



;i7) 



where i 



n— 1. 



Equation ([T7|) defines the transform [/ : {Btemp) -^ (-Bfrcq) appearing in the factorization 
Fg = t/ o (J (g) Ff/) o (p~^. What remains in order to obtain a concrete circuit computing the 
Fourier transform Fq, is an encoding of the bases for the transforms. With respect to the 
encoding 

ELeir') = k El^^C) = ^. 
the Fourier transform F^ for £,H defined by Equation (P) has the matrix representation Y^j. 
Let therefore the encoding of ^Btime be given by E'timelcV*^) = 2nj + k, and the encoding of 
-Btemp be given by £'temp(c-' ® C) = 2nj + i. With respect to this encoding, the transform 
(/ ® Fh) o ip~^ has the matrix representation I2 ®h F^^. 

Computing the U transform with respect to -Etemp is very simple 



\M) 



a i = or i = n 
ii i > n and i is odd 
otherwise. 



So, given a network computing the Fourier transform F = F2„ for the cyclic group of 
order 2n, a network computing the Fourier transform for the quaternionic group Qn can be 
constructed as follows. 



MSB 



LSB 



W-Z — MSB 



LSB 



Note that Z and W, both defined in Section |^, operate on distinct states and thus commute. 
If we, on the circuit given above, remove the Z gate, then we have a circuit that computes 
a Fourier transform for both the dihedral and the semidihedral group of order 4n. In the 
following, we show that this is no coincidence. 
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6.2 Metacyclic Groups 

In this section, we give a general quantum circuit for computing a Fourier transform for a class 
of metacyclic groups. A group is called metacyclic if it contains a cyclic normal subgroup H 
so that the quotient group G/H is also cyclic. Let G = {b^a^ '■0<j<q, 0<i< m} be a 
metacyclic group where 

b-^ab = a", 6" = a^ a"" = 1 

and with {m,r) = 1, m\s{r — 1), and q prime. Let d = {r ~ l,m). The group has a cyclic 
subgroup H generated by a of index q. Let T = {6^ : < j < g} be a left transversal for H 
in G, and write G = TH. Let -Etime : -Btime -^ '^qm be the encoding of Stime given by 



Va' 



nij + i. 



Theorem 14 The following network computes a Fourier transform for CG up to phase 
factors with respect to ii^time- Here uj = u^^. 



-^ 



m/d 



^[UJ] 



1-F — 



The phase factors involved in the theorem depend on the actual group structure. Before 
proving the theorem, we consider the representations of the group. The group G has qd one- 
dimensional representations, {p*-'}ij=o '■'"''" , each given by 

Let TZ^ be the complete set of representations of H given in Section |^. For every Q G TZ^ 
define the induced representation Q : G ^>- GLq(C) by 



iO^ 



OJ. 



\ri~ 



U)Z 



The group G has an if-adapted set of representations TZ consisting of the qd one-dimensional 
and (m — d)/q g-dimensional representations. The g-dimensional representations in TZ are 
all induced representations [|13 . 
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T®R^ 



The matrix coefficient p*-' G CG can be written as a linear sum of tlie basis-elements Btemp 

•H 
freq ' 

g&G 



^a^ 



keZq X£Zm 

^q \^qd 1^ ^m ^ tt ^ 

so, by definition of U : (Btemp) -^ (-Bfrcq) as given in Section ^ 

U-' : Vm,i ^ ^ E ^f (^S'^' ® &C'"/'',i,i)- (18) 

V y fcez. 

We refer to a matrix coefficient of an induced representation as an induced matrix coefficient. 
Any induced matrix coefficient (li G CG is non-zero on exactly one coset of H. For example, 
Cki — C31 is non-zero on the coset b^H. In general, the matrix coefficient Qi G CG can be 
written as a linear sum of the basis-elements ;Btemp as follows. 

Cki = Y.Cki(9)9 
geG 

teT h&H 

= T.cUt)T.Cuih)th 

teT heH 

= ai{h'~')Y.Qiih)h'-'h 

Here, (p = Cfc/(^'^ is some m-th root of unity. So 

U'^ : b^^,,l I — > 6''"' (g) b^.ri^i^y (19) 

that is, &^»,fe,i G i3freq is mapped by U~^ to one of the basis-elements i^temp up to a phase- 
factor. To find an expression for U instead of U~^, we need Lemma |l^ which easily follows 
from the following lemma for which a proof can be found, for example, in |TB, Lemma (47.8)]. 
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Lemma 15 The induced representation (^ is reducible if and only if there exists a j, 
1 < J < Q" — 1, such that ir^ = i (mod m). 



Lemma 16 Let C,li be any induced matrix coefficient. If C is irreducible, then ir^ is a not 
multiple of m/d. 

Proof To prove the contrapositive, suppose that ir' = {mod m/d). Since {r,m) = 1, 
i = (mod m/d), so id = (mod m,). Since d divides r — 1, we have ir = i (mod m,) and 
the statement follows from Lemma [T^. D 



Lemma 17 The transform U : (i3tcmp) -^ (-Bfrcq) is given by 

, \ (pbj if X is not a multiple of m,/d 

^ \ j=y,f Ejez, ^f &p-M,i if a: = tm/d. 

Here, is some m-th root of unity and bi G Bf^eq, both depending on the value of k and x. 

Proof Write ;Btemp as a disjoint union of two sets, ^B^emp ^^^ '^temp; where Bl^^^ = T ^ 
{&C^.i,i • 2; is a multiple of m./d}. Write similarly ;Bfreq as a disjoint union of two sets, B}^^^ 
and ;B|gq, where Bj^.^^ = {&pij,i,i : < i < d,0 < j < q}. We first show that 

{U-\BIJ) = (Bl^^) (20) 

by a simple counting argument. For each of the g(m — d) elements b^ij^i G Bf^.^^, we have 
that (^* G 7^ is irreducible. By Lemma ITO, ir' is not a multiple of m,/d, and therefore 



t/ ^{b^if^i) G (-Btemp)- Since iS^g^ and B"^^^^ has the same cardinality, and since U is unitary. 
Equation ( PP] ) follows. By Equation (|T^, the first case in the lemma follows . 

By the unitarity of U, we also have that 

(u-'iBL,)) = (^tU)- 

The action of U'^ on B}^^^ is given by Equation (|T8|), and the action of its inverse (that is, 
of U) on Bl is given by the second case in the lemma. D 



Let Ui : (Btemp) -^ {Bfreq) dcuotc the unitary transform which acts on Bl^^^ as U, and which 
^temp is given by 6^"' O 6^..^ ^^^ 1 — > b^.^^^i. 



on S^emp is given hy b'' ' ® b^i^i , , 1 — >• ft^i ^ ;. Since we are only interested in a quantum 
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network that computes a Fourier transform for CG up to phase factors, by Lemma O, it 



suffices to implement Ui instead of U. In conclusion, we have shown that the transform 

Fg = f/i o (/ ® F„) o (^-1 

is the Fourier transform for CG on 7?. up to phase factors. Here, Fm = Fh is the Fourier 
transform for Cif defined in Section ^ 

We now consider the implementation of Fq- The encoding -Etime : -Btime -^ '^qm is given 
above. Let -Etemp : -Btemp -^ '^qm be given by V ® C i— > mj + i. With respect to -Etime 
and -Etemp, the transform (/ ® F^) o ip"^ is implemented by I, ®r F^. Let the encoding 
-Efrcq of the matrix coefficients arising from the one-dimensional representations be given by 

p*-' t-^ jm + im/d. 

With respect to -Etomp and -Efreq, the transform Ui can be represented by 

I \km + im/d + x) if 1 < s < m/d 

\km + im d + x) \ — >< ^^^ ^ _ , ;. , . , • / , , \ .. n 

[Wgd ;7^ Ejez, ^f |J"^ + ^"^/f^ + x) if a; = 0. 

Here, /c G Zg, i G Z^, and x G Z^/t^. 

Written as a generalized Kronecker product, this is 

Fg ®iJ Id) X $gd(c^grf), Igd, . . . , Igd) ®i? Im/d- 

Thus, with respect to -Etime and -Efreq, -Fb is computed up to phase factors by a quantum 
circuit implementing 



F^ = (((Fg ®R Id) X ^qd{^ld)^ \d, • • • . V) ®/? Im/rf) ^ (iq 



JR ^ m 



Theorem 14 follows. 



7 Fourier transforms related to error-correction 

In this section, we give a quantum circuit for computing a Fourier transform for a certain 
subgroup En of the orthogonal group 0(2") = {A G GL2n(C) : AA* = I }. The group En was 
used independently by Gottesman |]18[ and Calderbank et. al. to give a group theoretical 
framework for studying quantum error- correcting codes. 

For alH = 1, . . . , n, define 

Xj = l2i-l ®R X ®R l2n-i 
Zj = l2'-i ®_R Z ®ij l2n-i 
Yj = l2«-l ®R Y (g)iJ l2n-i, 
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where X, Z, and Y are given as in Section ^. The group E^ is the group generated by 
these 3n unitary matrices. Its order is 2 ■ 4". Every element squares to either I or —I, and 
two elements either commute or anti-commute. When n = 0, En = {[±1]} is a cyclic group 
of order two, and if n = 1, En is isomorphic to D^. For larger n, En is isomorphic to D^/Kn 
where Kn is a normal subgroup isomorphic to Zg"^. Given a,c& Zg, a = (ai, . . . , an) and 
c = (ci, . . . , c„), let X(a) and Z(c), respectively, denote the elements IliLi X"* and HiLi ^iS 
respectively. Then every element g of En can be written uniquely in the form 

g = (-I)^X(a)Z(c) (21) 

where A G Z2, and a, c G Zg. We denote g by the 3-tuple (A,a,c). By rewriting Equa- 
tion (0), (7 can be written as a right Kronecker product 

g = (A,a,c) = ((-I2)'X'^^Z^^) 0^ (X'^^Z^^) ®« ■ ■ ■ »« (x""Z^"). (22) 

For n > 1, let i7 ^ En be the subgroup {(A,a, c) G -E„ : a„ = c„ = 0} of index 4, and 
identify E„_i with H in E„. Write E„ = TE„_i where T = {X^"Z^" : a„, c„ G Z2} is a left 
transversal for En-i in i^^- The group En has a complete set 7l(n) of 1 + 2^" inequivalent, 
irreducible and unitary representations, all but one of dimension one (except for n = where 
both representations, denoted ^^^p and ^^^a, are one-dimensional). The 2^" one-dimensional 
representations {^'^^p^^}x,zez" are given by 

Hp-^(^) = (")p"^((A,a,c)) = (-1)--+-^ 

The last representation, '^"V, has dimension 2" and is the group itself. From Equation (P^), 
we have the following recursive expression for the {kkn, lln)-th entry of the element g = 
(A, aan, ccn) G En, a, c,k,l e 1^2~^ , 

("V,,„,„((A,aa„,cc„)) = (-1)'"^" 5,„,„ ("-^)a,,((A, a, c)) 

where dn = kn®ln & Z2. Hence, 7^(„) is i?n-i^adapted relative to IZ^n-i)- 

We use the concept of adapted representations to find a Fourier transform for 'CEn- Let 
the bases ;Btime, -Bfreq, -B/fjng, S|^eq5 ^ud i^temp = T ® B^^^ be given as in Section ^. Let 
V? : (T ® i3/fjjjg) — i> (;Btime) dcuotc the natural isomorphism defined in Section E. 
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The matrix coefficients of 7l(n) can be written as linear sums of the basis-elements Btemp 



(n) XXnZZn 



E E ^"V""^^"((A,a,c))(A,a,c) 



an£Z2 c„eZ2 



'^kk„iu = Y. Y. ^^kk^iiAi^^ «' c)) (A, a, c) 

AeZ2 a,ceZ" 



E (-1)"^'" f E E ^""'kK(A, a, c)) (A, ad^, ccr.) 
cneZ2 V AeZ2 a,ceZ2~^ 

E(-l)^"'"^(X^"Z-®("-i)a,J 



where x, z,k,l E Z2 and dn = kn® In & '^2- Hence 



an&2 C„eZ2 






Cn6Z2 



(23) 



where fc„ = a„ © /„ G Z2. 

Equation (|23| ) seems to have the form of two W transforms for the one-dimensional represen- 
tations p, and a single W transform for the a representation. With respect to an appropriate 
encoding, this is indeed the case. Choose the encoding E^""^ : En — ^ Zg""^^, n >0, 

41((A,e,e)) = A 

-^time((A, aOn, CCn)) = E^^^^ ((A, a, c))anCn 
-C/templ-^n ^n '^ ^kl) — %req I ^kU'^nCn 

<(^°w) = 1 

/T-W /(n) xXnZ^nN _ j:p(n-l)/(n-l) xz\ 
-^freql T J ~ "'^freq I P J-^n^n 

4:i(^"K.„z/J = Et^^\^-~%,,)aX (where a^ = A;„ © /„). 

On the right hand side of the expressions, think of the images of the encoding as binary 
strings with standard string concatenation. 
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With respect to this encoding, the transform U : (-Btemp) -^ {^ircq), for n > 1, can be 
represented by 



|Asa„c„) 



I E.„ez. E.„ez.(-l)""^"+^"^"|Asx„^„) if A = 



where A G Z2, s G Zg" and a„, c„ G Z2. As a generahzed Kronecker product, this reads 



(I2 ®R (l22n-2 ®^ W, l22n-l)) ®R W. 



(24) 



For n > 1, let E be a quantum circuit computing the Fourier transform for CEn-i on 7^(n„i) 
with respect to the above encoding. Then, by Equation (^), the foUowing network computes 
the Fourier transform for CEn on ??.(„), also with respect to the above encoding. 



A — [ 



E 



W- 



W 



For n = 0, the one-bit network consisting only of the W transform computes the Fourier 
transform. Thus, expanding this recursively defined network given above, we have 

Theorem 18 The following network computes a Fourier transform for CEn- 



A -\W 
tti 

^^- 

an— 



W 



W- 



W 



w 



8 Conclusion 

The problem of finding efficient quantum algorithms computing a given unitary transform 
can be formulated as a purely matrix factorization problem. Let W be a set of basic unitary 
matrices. Given a unitary matrix U of dimension (nxn), can U be factorized into a product 
of basic unitary matrices such that the number of components in this product is polynomial 
bounded in log(n)? Previously, the only operations considered allowed in this product have 
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been the basic binary matrix operations: multiplication and standard Kronecker product. 
In this paper, we have shown that allowing a generalization of the latter, efficient networks 
can still be obtained. 

This generalized operation has several advantages. First of all, it gives a new tool when 
searching for factorizations of unitary matrices. The two new quantum networks given in 
Section ^]I| implementing the wavelet transforms were found this way. Secondly, it gives 
a nice compact mathematical description of more complex transforms. Thirdly, it directly 
gives quantum networks for computing unitary transforms which already were known to be 
expressible by generalized Kronecker products. This is for example the case for the Fourier 
transforms for the finite Abelian groups. 

In this paper, we have also discussed the issue of computing Fourier transforms for finite non- 
Abelian groups. We have given a definition of such computations on quantum computers, 
and especially we have given a slightly relaxed definition where we only compute a Fourier 
transform up to phase factors. Using this latter definition, we have devised a quantum 
network computing a Fourier transform for a class of meta-cyclic groups — even without 
completely knowing the group structure. This relaxed definition is in particular useful if the 



computation are to be followed by a measurement |19], as for example in the algorithms of 



Deutsch and Jozsa [|r^, Simon [^], Shor |2^, and Boneh and Lipton 0]. 

We have also given a simple quantum circuit computing a Fourier transform for a certain 
group |TB|, 0] used in quantum error-correcting. Together with Beals' proposal of a quantum 



network for the symmetric group [g], this emphasizes a challenging question which has only 
been partly discussed in this paper. Namely, which applications are there for these new 
transforms? Clearly, one can define quantum versions of the classical applications, but are 
there any other applications? For example, a crucial insight in Shor's algorithm |^9|] was the 



possibility of using the quantum version of the discrete Fourier transform to find the index 
of an unknown subgroup in a cyclic group. No efficient classical counterpart of this idea is 
known. Is this phenomena present for non- Abelian groups, too? 
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